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Propagation of light through media with a complex refractive index in which 
gain and loss are engineered to be PT symmetric has many remarkable features. 
In particular the usual unitarity relations are not satisfied, so that the reflection 
coefhcients can be greater than one, and in general are not the same for left or right 
incidence. Within the class of optical potentials of the form v{x) = vi cos(2/32;) + 
iv2sm{2/3x) the case V2 = vi is of particular interest, as it lies on the boundary of 
PT-symmetry breaking. It has been shown in a recent paper by Lin et al. that in 
this case one has the property of "unidirectional invisibility" , while for propagation 
^ . in the other direction there is a greatly enhanced reflection coefficient proportional 

Q I to L^, where L is the length of the medium in the direction of propagation. 

I For this potential we show how analytic expressions can be obtained for the various 

' transmission and reflection coefficients, which are expressed in a very succinct form 

in terms of modified Bessel functions. While our numerical results agree very well 
. with those of Lin et al. we find that the invisibility is not quite exact, in amplitude 

or phase. As a test of our formulas we show that they identically satisfy a modified 
version of unitarity appropriate for PT-symmetric potentials. We also examine how 
^ ' the enhanced transmission comes about for a wave-packet, as opposed to a plane 

Q ■ wave, finding that the enhancement now arises through an increase, of 0{L), in the 

. pulse length, rather than the amplitude. 
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PACS numbers: 42.25.Bs, 02.30.Gp, 11.30.Er, 42.82.Et 



INTRODUCTION - PT SYMMETRY AND OPTICS 



^ ■ The study of quantum mechanical Hamiltonians that are PT-symmetric but not 

Hermitian[l|-j^ has recently found an unexpected apphcation in classical optics0]-[lH, due 
to the fact that in the paraxial approximation the equation of propagation of an electromag- 
netic wave in a medium is formally identical to the Schrodinger equation, but with different 
interpretations for the symbols appearing therein. The equation of propagation takes the 
form 

where ip{x, z) represents the envelope function of the amplitude of the electric field, 2; is a 
scaled propagation distance, and V{x) is the optical potential, proportional to the variation 
in the refractive index of the material through which the wave is passing. That is, V{x) oc 
v{x), where n = no{l + v{x)), with |f | <C 1. A complex v corresponds to a complex refractive 
index, whose imaginary part represents either loss or gain. In principle the loss and gain 
regions can be carefully configured so that v is PT symmetric, that is v*{x) = v{—x). 
The set-up to which this applies is illustrated in Fig. 1: 
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FIG. 1: Propagation in the z-direction through a material with refractive index n = nQ{l+v{x)). In 
the paraxial approximation the equation of propagation is formally identical to the time-dependent 
the Schrodinger equation. 



Most of the initial applications of PT-symmetry in optics were concerned with such a set- 
up, where the variation in n is in the transverse direction. However, recently attention has 
turned to a different set-up, where the variation in n is in the longitudinal direction, as 



shown in Fig. 2. In particular the paper by Lin et al. [16| falls into this category. In that 
case the paraxial equation is not needed, and the scalar Helmholtz equation itself is formally 
identical to the time-independent Schrodinger equation. That is. 
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FIG. 2: Propagation in the z-direction through a slab of medium of length L with refractive index 
n = no(l + v{z)). The relevant quantities are now the transmission and reflection coefficients. 



It is this problem that we wish to discuss in the present paper, for the particular case 
v{z) = |a^e^*''^, which was the most interesting case discussed by Lin et al.[l6|. In the 
next section we briefly review the approximation scheme used in that paper to calculate the 
transmission and reflection coefficient. Then in Section 3 we develop an analytic formulation 
of the problem, whose solution and derived transmission and reflection coefficients can be 
expressed as simple combinations of modifled Bessel functions. In Section 4 we present our 
numerical results. Our conclusions are presented in Section 5, which also includes an analysis 
of the scattering of a wave-packet, as opposed to a plane wave, and a demonstration that 
the modifled unitarity relation developed by Ge et al. is satisfled identically. 
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II. UNIDIRECTIONAL INVISIBILITY 

The paper of Ref. [l6| dealt with the general case v{z) = Vi cos(2/3z) + iv2 sm{2/3z), with 
V2 not necessarily equal to Vi. The calculation was performed in a particular approximation 
scheme, whereby E was written as 

E = £f{z)e'^' + £b{z)e-'('\ (4) 

and second derivatives of and Sb arising from the Helmholtz equation (|3]) were neglected. 
Moreover, only terms proportional to e^*'^^ were kept, while those proportional to e^^^^^ were 
dropped, on the grounds that they were more rapidly oscillating. With these approximations, 
and keeping only linear terms in the detuning S = /3 — k, one obtains the following coupled 
linear equations for Sf, Sb- 

where e = (— 1^2^, —^ivik, 5), whose length is A = \J~^— ^k'^i^f — 

These equations are readily integrated, to give the following transmission and reflection 
coefficients for left or right incidence: 

Tl = Tr = X^D 

Rl = [e{v,-V2)'sm\XL)]/i4D) 

Rr = [k\v, + V2ysm\XL)]/{4D), (6) 

where the denominator is D = \/ X? cos^(AL) + 5^ sin^(AL). 

Note that although the transmission coefficients are the same for left or right incidence 
the reflection coefficients differ markedly. This is because the medium is directional: it is 
not symmetric under the parity operation z — )■ — z, but only under the combination PT . 

The special case when 77-2 = n\ exhibits the most striking results. In that case A^ = 5^, 
so that D = X, and Eq. (|6]) becomes 

Tl = Tr = 1 
Rl = 

= evlsm\6L)/6 . (7) 

Thus for left incidence one appears to have perfect transmission and no reflection. It is 
even the case that within the above approximations the transmission phase is identically 
zero, which gave rise to the terminology "unidirectional invisibility" . Equally striking is the 
large enhancement of the reflection coefficient for right incidence, which has the form of the 
square of a sine function, whose maximum is proportional to L. That is, for right incidence 
the reflected power grows like L^. 



(5) 
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III. ANALYTIC SOLUTION 

Recall that the analogue Schrodinger equation is 



+ k\l + 2v{z))E = (8) 



where in the case V2 = fi it is convenient to write v{z) = |a^e^*''^ 
Changing variables to y = (A;a//3)e^*'^^, the equation becomes 

This is the modified Bessel equation, with solution E = CI^{y) + DKu^y), where v = k/ (i. 
In the language of quantum mechanics {ip = E), 

,lj{z)=ClM + DK,{y) (10) 

This has to be matched onto ij) = A^e^^^ + B±e~^''^ at z = ±L/2 

A. Left Incidence 

Referring to Fig. 2, where the transmitted amplitude is normalized to 1 for convenience, 
we have the initial conditions 

^(z) = 1 =CUy+) + DK,{y+) 

^\z) = tk = {CIl{y+) + DK'M) X {t(3y+) (11) 

at z = L/2, where where y± = ua e^'^P^/'^, The solution for C and D is 

C = y+K^+i{y+) 
D = y+Iu+i{y+) , 

so that 

^{z) = y+ [K+^{y^)Uy) + h^,{y^)K,{y)] . (12) 
The first of Eqs. (ITTi) is satisfied by virtue of the Wronskian identity [3] 

K,+i{y)h{y) + h^,{y)K,{y) = 1/y . (13) 



Similarly 



ij'iz) = y^{iPy)[K,+,{yMy) + h+i{y+)K{y)] 

= [Ku+i{y+)lu+i{y) - lu+i{y+)K^+x{y)] 

+^ [K,+,iy+)lM + h+i{y+)Kiy)] }, (14) 
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which satisfies the second of Eqs. (fTTi) by virtue of the same identity (recall that k = 1^(3). 
At z = —L/2 we have to match ip with A^e^^^ + B^e'^^^ . The general formulas are 
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which give, after some algebra, 

Bl = (^«V)e-^'=^[-J^.+i(i/+)/.+i(l/_) . (15) 

B. Right Incidence 

The set-up is shown in Fig. 3, with the transmitted amplitude again normalized to 1. 

L 
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FIG. 3: Set-up for propagation from the right 



The initial conditions at z = —L/2 are 

^(z) = 1 =CI,{y^) + DK,{y.) 

^'(z) = -rk = {CIliy^) + DKliy.)) x {tf3y.) . (16) 
The solution for C and D is 

C = y-K,_,{y_) 
D = y_I^_i{y_), 

so that 

^{z) = y. [K,.,{y.)lM + h.,{y.)K,{y)] . (17) 

Similarly 

i^'{z) = y4t(3y)[K,^i{y^)Iliy) + U^iiy^)K{y)] 

= y.m)\ [K-i{y-)h-i{y) - h-i{y-)K,-M] 
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^ [K,^,{y^)lM + Iu^i{y^)KM] }• (18) 



Again, the initial conditions f ll6p are satisfied by virtue of Eq. f ll3p . 

At z = L/2 we liave to matcli ip witli A^je"*'"'^ + B^e^^^. Tlie general formulas are 
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giving 



Ar = {\a^u)e'^'^[ 4_i(y_)JC+i(z/+)-/C-i(l/-)/.+i(y+)] 

5^ = (^«V)e-^'=^[-J._i(y_)K,_i(y+) + K,_i(y_)J,_i(y+)] . (19) 

C. Analytic Continuation 

Equations f|T5|) and f|T9l) constitute our main results, but before going on to the numerics 
a word of caution about the evaluation of the modified Bessel functions is in order. Recall 
that the argument is ?/ = [ka / [5)e^'^^^ . As z goes from —L/2 to L/2 this encircles the origin 
many times, crossing the cut on the negative real axis. Thus we need to know how to 



continue onto subsequent sheets, which is achieved by using the continuation formulas [18 

Kiye^n = e-^^-KM - ^^^^^/.(y) • (20) 

sin(7rz/j 

The resulting functions are smooth functions of with no discontinuities, as illustrated in 
the following diagram for lo^iv)- 
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IV. NUMERICAL RESULTS 



The parameters we use in this section are those of Ref. [16|, namely uq = 1, rii = 0.001, 
L = 12.5%, (3 = 100. 



A. Left incidence 



For left incidence, we find that is very close to 1, over a wide range of 5, so that the 
transmission coefficient T = l/|ylj;^p is very nearly 1, in agreement with Ref. [l6|. On closer 
inspection, however, there is a very small deviation of RelAi) from 1 and /m(y4i) from 
near 5 = 0, as shown in Fig. 5. 
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FIG. 5: (color on line) Real and imaginary parts of Al as functions of S. 



Correspondingly the transmission is not quite perfect, in amplitude or phase, as shown in 
Fig. 6. 
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FIG. 6: Transmission coefficient and phase of the transmission amplitude as functions of 5. 



The reflection coefficient Rl = \Bl\'^ /\Al\'^ is indeed very small, but varies rapidly on the 
scale of 10~^, as shown in Fig. 7. 
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FIG. 7: Reflection coefficient Rl function of 5. 



B. Right incidence 



For right incidence the transmission properties are identical, since A^i = A^, but the 
reflection coefficient differs markedly, as expected from Eq. ((Tj). In Fig. 8 we show the real 
and imaginary parts of Bf> and the corresponding reflection coefficient Rr = [-BrP/IArP 
The latter is almost indistinguishable from the corresponding figure in Ref. 
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Right panel: reflection coefficient Rr as a function of 5. 



V. DISCUSSION 



Our main results are the simple analytic expressions, Eqs. f fTSj) and f|T9|) . for the trans- 
mission and reflection amplitudes in the case V2 = vi. It transpires that the approximate 
expressions given in Ref. [iGj are surprisingly good, but that the transmission amplitude 
A^iR and the reflection amplitude Bl exhibit interesting oscillatory behaviour on a small 
scale. The large reflection coefficient Rr., which is perhaps the most striking effect of the PT- 
symmetric refractive index, is essentially indistinguishable from the simple form of Eq. 

A test of our expressions is afforded by the modified unitarity relations recently found by 
Ge et al. for PT-symmetric potentials. Instead of the usual unitary relation T + R = 1, 
applicable to a Hermitian potential, one now has the equation 



T-1 = ±^/RlRr 



(21) 
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In the present case, as can be seen from Fig. (6), the transmission coefficient T is greater/less 
than 1 for S positive /negative. We can write each side in terms of the amphtudes A, B^ ^ 
themselves without modulus signs by noting from Eqs. (IT^ and (IT^ that the complex 
conjugates of Ae~^^^ and Br^lc'"^^^ are obtained by interchanging the arguments y±. This 
shows that Brc^'^^ and Bie^^^ are pure imaginary, so that \BrBl\ = ±(i?/je*'^^)(i?Le*^-^)* 
depending on whether Br and Bl are in or out of phase. After some algebra it can be seen 
that the modified unitarity relation is satisfied identically by virtue of the Wronskian-type 
identity 

- = ^, (22) 

y 



which is not found in Ref. [18|, but is readily proved. 

In the more general case V2 < vi, an analytic solution is still possible, in terms of Mathieu 
functions with complex argument, using the similarity transformation jl9| that transforms 
the potential t"! cos(2/32;) + iv2sin(2Pz) into the Hermitian potential y/{vl — v^) cos{2Pz). 
However, from the mathematical point of view the resulting expressions do not simplify in 
the manner of Eqs. (ITS!) and (fTOl) . because of the lack of recursion relations for the Mathieu 
functions, and from the physical point of view the results are not as dramatic. 

Finally it is interesting to investigate how the enhancement of the reflection coefficient 
is affected when the input is not a continuous plane wave, but rather a wave packet of finite 
width in z or t. The question is of interest because on the one hand it has been shown [20[ on 
general grounds that such an enhancement is a general feature of PT-symmetric potentials 
at the symmetry-breaking point, while on the other hand, in the diffraction set-up of Fig. 1 



it has been shown [12|, |21| that for a wave-packet input the maximum amplitude becomes 
saturated, and the enhanced output is obtained instead from a spreading of the beam. 
Including the time dependence, the field ipki^, t) for z > L/2 in the case of right incidence 

is 



where for n very close to 1 we can take u = k in natural units. In Eq. fl23l) we may take 
Ap^^I and Br^ \ia^ke'^^ sin{5L) /5 to a very good degree of approximation, in accordance 
with Eq. d?]), so that the right-going component is 

ij^{z,t) = l^a2g./3L^^M-V2-t)!^. (24) 
2 



Then, if we modify the incoming component of the field by folding with a Gaussian 
shape (l/(w;0F)e-('=-^)'/'"', the left- going component t) becomes simply a Gaussian 

proportional to e"**" {z+L/2+t) ^ ^j^ijg (recall that k = (3 — 5) the right-going component 
becomes 



oo 



in which we can neglect 5 in the factor (/3 — 5) occurring in the integrand, since the large 
enhancement comes from the /3 term alone. 
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By using the integral representation | J^j^ dne^^^ for sin 6L/6 and interchanging the orders 
of integration we can write 



d5 e-^'/"''+*'^« 



sin 5L 

6 



^vr [erf + o) + erf - 



(26) 



in which ^ is to be identified with z — L/2 — t. For a wide range of w the two error functions 
behave hke sgn functions, and the combination of the two gives an almost constant plateau 
for 1^1 < L, while essentially vanishing outside that range, as is illustrated in Fig. (9) for 
w = 1 and L = 10. Thus for a wave-packet input the enhanced reflection comes not from an 
increase in the overall amplitude, but instead from a lengthening of the pulse, of order L. 
An analogous effect was noted in the diffraction situation by Longhi 12[, and the connection 
to the properties of error functions pointed out in Ref. 
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FIG. 9: Generic shapes of input and output pulses (left and right panels respectively). 
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